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Abstract. This Chapter outlines the replica approach in Random Matrix Theory. 
Both fcrmionic and bosonic versions of the replica limit arc introduced and its trickery 
is discussed. A brief overview of early heuristic treatments of zero-dimensional replica 
field theories is given to advocate an exact approach to replicas. The latter is pre- 
sented in two elaborations: by viewing the (3 = 2 replica partition function as the Toda 
Lattice and by embedding the replica partition function into a more general theory 
of T functions. 
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1. Introduction 

1.1. Resolvent as a field integral 

In physics of disorder, all observables depend in highly nonlinear fashion on a stochastic 
Hamiltonian hereby making a nonperturbative calculation of their ensemble averages 
very difficult. To determine the average quantities in an interactionless system, one has 
to know the spectral statistical properties of a single particle Hamiltonian contained 
in the mean product of resolvents G{z) = tr — iK)~^ defined for a generic, complex 
valued argument 2; G C \ M. Each of the resolvents can exactly be represented as a ratio 
of two functional integrals running over an auxiliary vector field i/j which may consist 
of either commuting (bosonic, 1^ = 3) or anticommuting (fcrmionic, = X) entries. 

In the random matrix theory (RMT) limit, when a system Hamiltonian (or a 
scattering matrix) is modelled by an N x N random matrix 3-C of prescribed symmetry, 
the resolvent G{z) equals 

G{z) = ir]S,J ipeipee''^^^^''"''''!'^ (^J D['0,iA] e^"^^'^^"''^''^)^ . (1.1) 

Here, S^h: = ''/'(-^ — IK)i/? = 'ipj {zSjk — "Kjk) ipk (summation over repeated Latin indices 
is assumed throughout this Section); vector i/j is defined by V = {i^i: ■ ■ ■ :'4'n)'^ , 
whilst = (■01, • • • , ■0jv) is its proper conjugate. The parameter 77 accounts for 
the nature of both vectors: it is set to +1 for fermions {ip — x) ^^^^ to —1 for 
bosons (i/j = s). In the latter case, convergence of field integrals is ensured by the 
regularizer = sgn(3m.z). The notation !D['0,'0] stands for the integration measure f 

a5[V',^] = (27r)-^nf=i^V^r 

Equation (1.1) may conveniently be viewed as a consequence of the identity 

G(z)=r)-^\ogdet%z-3£), (1.2) 

combined with the field integral representations of the determinant (?] = +1, -0 = x) 
and/or its inverse (77 = —1, if) = s): 

det''(^- JC) = {ivSzf J e^^^^«("'*'^\ (1.3) 

Although exact, both bosonic and fermionic versions of Eq. (1.1) are a bit too 
inconvenient for a nonperturbative ensemble averaging due to the awkward random 
denominator. 

To get rid of it, several field theoretic frameworks have been devised by theoretical 
physicists: the replica trick [Weg79, SchSO, Efc80], the supersymmetry method 
[Efe82a, Efe82b], and the dynamic (Keldysh) approach [Hor90, Kam99a]. Leaving aside 

t Integrals over anticommuting (Grassmann) variables x and x are normalised according to 

jdxx = j rfxX=(27r)V2. 
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the supersymmetry and the Keldysh techniques this Chapter aims to provide an 
elementary introduction to notoriously known replica approach whose legitimacy has 
been a point of controversy [Ver85b, Kam99b, Zir99] for over two decades. Recently 
discovered integrability [Kan02, Kan05, Osi07] of zero-dimensional replica field theories 
will be a dominant motif of this contribution. 



1.2. How replicas arise and why they are tricky 

Both fermionic and bosonic replicas are based on the identity § 

logX=lim^— ^ (1.4) 

n^O n 

which can be very useful in evaluating the average of a logarithm (log X) of the random 
variable X. Indeed, identifying X with det ''{z — IK) in Eq. (1.2), and further combining 
Eq. (1.4) with the field integral representation Eq. (1.3) of the (inverse) determinant, 
we formally relate the resolvent 

d 

a — 

'0 n 

to the 'partition function' 



n p 

2;W(^) = n / '^['^^"'\ '^^°'^] exp [is, '^^'"\z - 5<:)V'^''^] (1.6) 

a=l'' 

of n copies, or replicas, of the initial random system. The pair of formulae Eqs. (1.5) and 
(1.6), known as the replica trick ||, achieve our goal of removing a random denominator 
from Eq. (1.1) and hint that a nonperturbative calculation of the mean product of 
resolvents may become feasible. 

In order to keep the discussion concrete and set up notation, we further assume that 
the matrix Hamiltonian is drawn from the paradigmatic Gaussian Unitary Ensemble 
(GUE) associated with the probabihty measure ^ [Meh04] 

N N 

j=i i<j<k 

Then, the resolvent g{z) = {G{z)) averaged with respect to the GUE probabihty measure 
should be furnished by the limiting procedure 

I An introductory exposition of the supersymmetry method can be foimd in Chapter 7 of this Handbook 
as well as in the earlier review papers [Efe83, Ver85a] and the monograph [Efe97]. For a review of the 
Keldysh approach, the reader is referred to [Kam09]; see also the paper [AltOO] where the Keldysh 
technique is discussed in the RMT context. 

§ Mostly known in condensed matter physics community since the paper by Edwards and Anderson 
[Edw75] on spin glasses (see also [Eme75]), the recipe of calculating the average of a logarithm based 
on Eq. (1.4) dates at least as far back as 1934, see the book [Har34] by Hardy, Littlewood and Polya. 

II The alert reader might already detect some trickery behind Eqs. (1.5) and (1.6). Their mathematical 
status will be clarified in Section 1.2.3. 

^ See also Chapter 4 of this Handbook which outlines the method of orthogonal polynomials for unitary 
invariant random matrix models. 
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where the (average) replica partition function equals 

(4^n^)) = n / exp(-i,75.ztro- + |trfr^). (1.9) 

a=l 

Here, the Hermitean matrix cr — a"^ acting in the rephca space is defined as 

a^p^{i^j®i^j)^p^^fi^f\ Q;,/3e(l,...,n). (1.10) 

Notice that ensemble averaging of Zn\z) [Eq. (1.6)] has induced an efi^ective 
interaction between n replicated random systems as described by the interaction term 
tr<T^ = tr [('0J ® V^j)^] in the action. To make further progress in evaluating {Zn\z))^ 
the interaction term may routinely be decoupled by means of the Hubbard-Stratonovich 
transformation [Efe97, Ver85a]. 

1.2.1. Fermionic replicas 

In case of fermionic fields (V' = X) ^ = +1); the x'^ interaction can be decoupled 
via n X n Hermitean matrix field Q„ as 

expl+hraA = tt-^'/^ / 2)Q„ exp [-tr (O^ + orO„)] . (1.11) 
Performing the integration over fermionic fields in Eq. (1.9) with the help of Eq. (1.3), 
n / expf-is^ztrcT + h^TtrA 

a=l V / 

= 7r-"'/2(-is,)"^ / 2)Q„ e-*'^^" det^(^ - isM, (1.12) 
we arrive at the fermionic replica limit 

^(z) = lim 1^(2;W(^)) (1.13) 

n^O n oz 

that relates the average resolvent g{z) to the fermionic replica partition function 

{Z^+\z)) = [ 2)Q„ e-*'-^" det^ {z - is,Q„) . (1.14) 
Before discussing this result, let us turn to the bosonic version of the rephca trick. 

1.2.2. Bosonic replicas 

In case of bosonic fields = s, rj = —1), decoupling of the interaction can be 
carried out in a similar fashion. Making use of yet another variant of the Hubbard- 
Stratonovich transformation 



exp (~tr crA = tt-"'/" / l>0„exp [-tr (O^ + i(TQ„)] , (1. 

\ 4 / JQl=Qn 



15) 
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and integrating out bosonic fields in Eq. (1.9) with the help of Eq. (1.3), 




a=l 



= 7r-"'/2(+w,)"^ / DQ„ e-^-^O" det-^(z - 5,Q„), 



(1.16) 



jQi=Qn 

we express the average resolvent 



(1.17) 



through the bosonic replica partition function 



{Zi-\z)) = / D0„ e-*--^" det-^ {z - 5,0„) . 



(1.18) 



jQi=Q, 



1.2.3. Subtleties of the replica limit 

Seemingly innocent at first glance, both fermionic [Eqs. (1.13) and (1.14)] and bosonic 
[Eqs. (1.17) and (1-18)] replica prescriptions appear to be counterintuitive and rising 
fundamental mathematical questions [Par03]. Indeed, due to a particular integration 
measure which makes no sense for n other than positive integers {n e Z+), the matrix- 
integral representation of average replica partition functions 



cannot directly be used to implement the replica limit Eqs. (1.8) as the latter is 
determined by the behaviour of {Zn\z)) in a close vicinity oi n — 0. This mismatch 
between the 'available' {n e Z"*") and the 'needed' (n G M) is at the heart of the trickery 
the replica field theories have often been charged [Ver85b, Zir99]. 

The canonical way to bridge this gap is to determine the average replica partition 
functions {Zn\z)) for n G and then attempt to analytically continue them to n G M, 
in general, and to a vicinity of n = 0, in particular. This is a nontrivial task for two 
major reasons: (i) The analytic continuation of the replica partition function away from 
n integers should not necessarily be unique (ii) To retain control over the analytic 
continuation, the latter must rest on an exact calculation of the average replica partition 
function for n G Z+. Early approaches to rephca field theories seem to underestimate 
these two points bringing a number of pathological results even in the RMT setting (see 
Section 2). 

2. Early studies: Heuristic approach to replicas 

Exact evaluation of replica partition functions, whilst welcomed, is quite a challenge. 
At the same time, their approximate calculation is often feasible in a certain region of 

^ For entire functions, the uniqueness is guaranteed by a boundedness property as formulated by 
Carlson's theorem [Tit32]. 




(1.19) 
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parameter space where a saddle point procedure can be justified. In doing so, one is 
naturally led to the 'replica symmetric ' and 'replica asymmetric ' saddle point manifolds 
as discussed below. 

2.1. Density of eigenvalues in the GUE 

In the RMT context, a saddle point evaluation of the replica partition function Eq. (1.19) 
makes sense if the dimension N of the random matrix IK is large enough. For not 
too large replica parameter n G Z+ (in particular, n should not scale with N), the 
dominating contribution to {Zn\z)) is expected to come from the configurations Q^^''^ 
determined by the saddle point equation 

— -trfe - r/iV log(^ - Vi^^\,QA = 0. (2.1) 

Its solutions form 2" saddle point manifolds 

with taking on the values ±1 independently of each other. Here, 

e''^ ^z, + i^Jl-zl (2.3) 

where Zs stands for the scaled energy Zs — zfDedge with Dedge = V2N being the endpoint 
of the spectrum support. (Hence, the spectrum bulk is situated within the segment 
|9^e^s| < 1-) 

2.1.1. Bosonic replicas 

Out of the plethora of saddles Eq. (2.2), only the distinguished replica symmetric man- 
ifold 

^t^^ ^^.Me^'®tn (2.4) 

sym V 

contributes the bosonic rephca partition function Eq. (1.18). This is so because Eq. (2.4) 
is the only saddle [Zir99] (i) reachable by continuous deformation of the integration 
contour in Eq. (1.18) without crossing the hypcrsm^face defined by the singularities of 
det^^(2; — SzQ„) and (ii) compatible with analyticity of the average resolvent Eq. (1.17) 
at infinity. In the leading order in the large parameter iV, the bosonic replica partition 
function is then approximated by 

(2;H(. = e - iO)) ^ (^^y (2sin ey^'/^e^p (e^^^ - - 4 " £ 

(2.5) 

By derivation, Eq. (2.5) holds for n G Z"*". To retrieve the average density of 
eigenlevels 

^(e) = -Jm^(e-iO) (2.6) 

TT 
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through the rephca hmit Eq. (1.17), one should analytically continue {Zn {z)) away 
from n e Z"*". To be on the safe side, such an analytic continuation must rest on an 
exact integer-n result for {Z^\z)). The latter is sadly unavailable. Not being spoilt 
for choice, one could try to analytically continue the bosonic replica partition function 
to the domain n e R"^ taking the approximate result Eq. (2.5) as a starting point and 
merely assuming it to hold, as it stands, in the right vicinity of n = 0. Then, the replica 
limit Eq. (1.17) can be taken to yield the Wigner semicircle [Edw80] 

^(es) = -yr^, |es|<l. (2.7) 

71 

Here, eg = e/Dgdge- The 1/N correction to Eq. (2.7) can be obtained in a similar fashion 
and is known to vanish within the replica symmetric ansatz. 

Similarly, replica symmetric saddle point calculations performed for the Gaussian 
Orthogonal Ensemble (GOE) and the Gaussian Symplectic Ensemble (GSE) of random 
matrices yield [Edw76, Ver84, Dhe90, Ito97] 

esl < 1. (2.8) 

Here, (3 is the Dyson symmetry index [Meh04] taking the values /? = 1 for GOE, j3 = 2 
for GUE, and /? = 4 for GSE. 

Importantly, both leading and subleading inl/N terms in Eq. (2.8) for the average 
densities of eigenlevels contain no terms oscillating on the scale of the mean level spacing. 
Since the replica asymmetric saddles are inaccessible as explained on general grounds 
[Zir99] in (i) and (ii) below Eq. (2.4), we are led to conclude that the saddle point 
approach to bosonic replicas fails to reproduce truly non-perturbative features of the 
eigenlevel density. 



1 - 



- ^S{el - 1) 



2.1.2. Fermionic replicas 



Looking for some insight into a possible role played by replica asymmetric saddles, 
we turn to the approximate performance of fermionic replicas Eq. (1.14). In this case, 
all 2'^ saddle point manifolds 



Oi^P) = J diag(e-i^ . . . , e-"^), Ke = ±1, 



are accessible adding up, for n — 2m, to 



'N 



mN q=+m / at \ m^-q^ 



^mN cos 26 



(2 sin e) 



q=—m 



X 



V27ry 

exp ^iq (^N{2e - sin 26) + 2m(^e-^ 



Here, Vm,q denotes the volume of Grassmanian [Kam99b, Zir99] 



'^m,q = vol 



U(2m) 



U(m — q) X U(m + q) 



(27r) 



nr=i r(j) nr=7r(j) 



n?=ir(i) 



(2.9) 



(2.10) 



(2.11) 



2 Early studies: Heuristic approach to replicas 



9 



The summation index q in Eq. (2.10) counts (2m + 1) families of saddle point manifolds 
[Eq. (2.9)], the q-th family being represented by the configuration 

Qi^^Q] = -i5.y^diag(e^^ ® e'"' (g) 1^+,) (2.12) 

taken with the obvious combinatorial weight 
2m \ _ r(2m + 1) 



m + qj r{m + l-q)r{m+l + q)' 

Obtained in the large- limit, Eq. (2.10) holds for m e Z+ which do not scale with N. 

Aimed at deriving the density of eigenlevels through the replica limit, one should 
first to analytically continue Eq. (2.10) away from m e Z"*". Even though making an 
analytic continuation based on an approximate result is a dangerous ploy and is, with 
certainty, a mathematically questionable procedure, we embark on the proposal due to 
[Kam99b] who have spotted that the volume of Grassmanian [Eq. (2.11)] vanishes for 
all integers \q\ > m + 1. This observation makes it tempting to extend the summation 
over q in Eq. (2.10) to (minus and plus) infinities to end up with the following trial 
function for 'analytically continued' fermionic replica partition function 

/'AT\'^^ (?=+oo , .m^-q^ 

- (f ) e---^ y: ^rn, (^) (2sin er-^^' 

^ ^ q=—oo ^ ' 

X exp \iq (^N{29 - sin 29) + 2m (^^ - |) ) . (2.13) 

A close inspection of this result reveals that it is fiawed: (i) so-continued replica partition 
function (2,2^^(2;)) diverges + [Kam99b, Zir99] in the vicinity of m = -|-0, the region 
crucially important for retrieving the density of eigenlevels through the replica limit; 
(ii) due to the q ^-^ —q symmetry of the summand, {Z^2mi^)) must be real thus leaving 
no room for a non-zero density of states. Indeed, a formally derived small-m expansion 
of Eq. (2.13) 



(2;W(^))^l+m 



/V\ 1 

N ( cos 29 + log - j + ^^^p^ cos [Ar(2^ - sin 2^)] 



+ 0{m^) 

(2.14) 

considered together with Eq. (2.6) leads us to conclude that the replica limit [Eq. (1.13)] 
in the above elaboration fails to reproduce even the smooth part of the average density 
of eigenlevels yielding g^e) ~ 0. Notice that the derivation of Eq. (2.14) boldly ignores 
the divergence of the infinite series Eq. (2.13). 

This unphysical result is at odds with the work [Kam99b] where both the Wigner 
semicircle and the 1/N oscillating correction to it. 

Pie) ~ ^ [1 - ^ cos(Ar(2g-sin2g)) 1 

Q{e)-^V^ '4 4Ar sin^^ J' ^^'^^^ 

were reproduced out of fermionic replicas in almost the same elaboration. The only 
difference between the above treatment and that in [Kam99b] is that its authors used an 

The group volume Vm^q grows too fast with q for the series J2gZt'^{- ■ •) to converge. 
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alternative enumeration of saddle point manifolds Eq. (2.12) contributing the fermionic 
replica partition function. 

To meet the parameterisation of [Kam99b] , we introduce a new summation index, 
p — q + m, in Eq. (2.10). This amounts to counting saddle point manifolds starting with 
the 'replica symmetric' one, 



I'^Vl = diagte-" ® l2„-p, e"'" » 1,) (2.16) 



as p varies from to 2m, so that 

„mN cos 29 2m / , ^ \ p(2m-p) 



\.ip{2m—p) 



/ A[\™-"' pTuN cos 29 / \PK^rn-p) 

X exp [i(p-m) (^Ar(2^-sin2^) + 2m (^^- I)) . (2.17) 

For m e Z"*", Eqs. (2.10) and (2.17) are trivially identical. However, this is not the case 
for generic real valued m after Eq. (2.17) is 'analytically continued' by extending the 
summation over p to (minus and plus) infinities: 



/ „mN cos 29 P-+°° / j> r n p{2m-p) 

X exp [i(p-m) (^iV(2^-sin2^)+2m(^^-^)) . (2.18) 

Now, contrary to the previous result [Eq. (2.14)], the small-m expansion of the fermionic 
replica partition function appears to contain a mysterious imaginary component *, 



N { e^^^ - 2i0 + log -] + e^^(2^-^-2^) 



+ 0(m^). (2.19) 



4iV(sini 

This expansion coincides with the one in [Kam99b] and does reproduce, through the 
replica limit, the correct result [Eq. (2.15)] for the average density of eigenvalues in GUE, 
in both leading and subleading orders in 1/A^. As the latter term captures oscillatory 
behavior of the eigenlevel density, it was assumed in the literature [Kam99b, Kam99c, 
Yur99] that the replica asymmetric saddles may describe a nonperturbative sector of 
replica field theories. 

Let us stress that neither of the above two treatments of fermionic replicas [resulting 
in Eqs. (2.14) and (2.19), respectively] can be considered as mathematically satisfactory 
because they both rely on a nonexisting 'analytic continuation' of the replica partition 
function [Eq. (2.10)] to the vicinity of m = +0, as explained below Eq. (2.13). 



* It was argued in [Zir99] that the procedure of 'analytic continuation' that led to Eq. (2.18) and further 
to Eq. (2.19) favors so-called causal saddle points over their conjugate counterparts called acausal. Such 
a selectivity is eventually responsible for the correct result for the density of eigenlevels as discussed 
below Eq. (2.19). 
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2.2. Brief summary 

A brief tour d'horizon on the saddle point approach to rephca field theories indicates that 
the bosonic variation of the replica trick is restricted to the perturbative sector of the 
theory accounted for by the (only reachable) replica symmetric saddle point manifold. In 
the fermionic version of the replica trick, both replica symmetric and replica asymmetric 
saddle point manifolds contribute the replica partition function; however, its analytic 
continuation to a vicinity of n = is ill defined, both in terms of convergence and 
uniqueness. 

These drawbacks of the saddle point approach to replicas are not specific to the 
GUE description. Similar difficulties arise in replica studies of other random matrix 
ensembles [DalOl, Nis02] and in analysis of physical systems (notably one-dimensional 
impenetrable bosons) admitting effective RMT description [GanOl, Nis03, Gan04]. 

Are the problems surfaced in the above calculation indicative of internal difficulties 
of the replica method itself or should they be attributed to a particular computational 
framework? A little thought shows that the approximate evaluation of both bosonic and 
fermionic replica partition functions is the key point to blame for the inconsistencies 
encountered in the above elaboration of the replica trick. In such a situation, leaning 
towards exact calculational schemes in replica field theories is a natural move. 

3. Integrable theory of replicas 

In this section, we outline an alternative way of treating replica partition functions. 
A connection [Kan02, Kan05] between zero-dimensional repfica field theories and the 
theory of integrable hierarchies f is central to our formalism. 

3.1. Density of eigenvalues in the GUE revisited (easy way) 

For illustration purposes tJ, we choose the very same problem of calculating the average 
density of eigenlevels in the GUE specified by the probabifity measure Eq. (1.7). For 
the lack of space, only fermionic replicas will be considered; a bosonic replica treatment 
can be found in the tutorial paper [Osi09] . 

3.1.1. Replica partition function as Toda Lattice 

Our claim of exact solvability of the replica model Eq. (1.14) and the models of the 
same ilk rests on two observations. To make the first, we routinely reduce the average 

f An introductory exposition of integrability arising in the RMT context can be found in Chapter 10 
of this Handbook. 

)i Since the framework to be presented here solely rests on the symmetry underlying the matrix model, 
it can readily be adopted to other spectral statistics for random matrix ensembles falling into the same 
P = 2 Dyson's symmetry class [Meh04]. The reader interested in a general formulation of the integrable 
theory of replicas is referred directly to Section 3.2. 
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fermionic partition function Eq. (1.14) to the n-fold integral 



(3.1) 



after diagonalising the Hermitean matrix Q„ = UnXUl^ by unitary rotation G U(n); 
here A is a diagonal matrix A = diag(Ai, . . . , A„) composed of eigenvalues of Qn, and 
A„(A) is the Vandermonde determinant 



A„(A)=det[An = n(A,-Afe) 



(3.2) 



induced by the Jacobian of the transformation Q„ h-* {ILn, A). Further, making a proper 
shift of the integration variables and applying the Andreief-de Bruijn integration formula 
[And83, deB55] 



]^d//(A^) det[/fe(A^)]det[g'fe(A^)] = n! det 



dfi{X)fk{X)gi{X) 



which holds for benign integration measure diJ,{X), one derives: 

(4+)(^)) = exp [nz^] det k,e \ [ dA A^+'^+^exp (-A' - 2izX) 

The latter is equivalent to the remarkable representation 

(2;W(^))=exp[n^2] 
involving the Hankel determinant 



rW(z) = det arV(+)(z) 



(3.3) 



(3.4) 



(3.5) 



(3.6) 



with r['^\z) 



k/=0,-,n-l 

e~^'^ Hn{z) being related to the Hermite polynomial Hj^[z) (see also 
Chapter 4 of this Handbook). In the above equations, no care was taken of prefactors 
which tend to unity in the replica limit. 

Consequences of the Hankcl-dctcrminant-likc representation Eq. (3.5) of the 
fermionic replica partition function {Z[^\z)) are far reaching. As had first been shown 
by Darboux [Dar72] a century ago, any set of Hankel determinants meeting the 'initial 
condition' Tq'^\z) — 1 (which is indeed the case for Eq. (3.6) due to Eq. (3.5) and the 
normalisation {2!^\z)) — 1) satisfies the equation 

^2 /a X 2 



dz'^ 



—T^+Hz) 



r!C-\{z)TiX\{z), neZ+. (3.7) 



Equations (3.5) and (3.7) taken together with the known initial conditions for Tq'^\z) 
and Ti'^^ {z) establish a differential recursive hierarchy between nonperturbative fermionic 
replica partition functions {Zl^'\z)) with different replica indices n. This is an exact 
alternative to the approximate solution Eq. (2.10) presented in the previous section. 

Equation (3.7), known as the positive Toda Lattice equation [Tod67] in the theory 
of integrable hierarchies * [Mor94] , is the first indication of exact solvability of replica 



See also Chapter 10 of this Handbook. 
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field theories. Importantly, the emergence of the Toda Lattice hierarchy is eventually 
due to the (5 — 2 symmetry of the fermionic replica field theory encoded into the squared 
Vandermonde determinant in Eq. (3.1). 

3.1.2. From Toda Lattice to Painleve transcendent 

While important from conceptual point of view, the positive Toda Lattice equation 
for the fermionic replica partition function {Z^\z)), if taken alone, is not much helpful 
in performing the replica limit. 

Fortunately, here the second observation, borrowed from [ForOl], comes in. 
Miraculously, the same Toda Lattice equation governs the behaviour of so-called r- 
functions arising in the Hamiltonian formulation [Oka86] of the six Painleve equations ft 
[Cla03, Nou04], which are yet another fundamental object in the theory of nonlinear 
integrable systems. The Painleve equations contain the hierarchy (or replica) index n as 
a parameter. For this reason, they serve as a proper starting point [Kan02] for building 
a consistent analytic continuation of replica partition functions away from n integers. 

The aforementioned Painleve reduction [Oka86, ForOl] of the Toda Lattice equation 
Eq. (3.7) materialises in the exact representation [ForOl] 

(2;W(z)) = (2;W(0))exp (^£dtaMt)) (3.8) 

which holds as soon as n G Z+. It involves the fourth Painleve transcendent o"iv(t) 
satisfying the Painleve IV equation § in the Jimbo-Miwa-Okamoto form [Jim81, Oka86] 

{a'^^f - A{ta{^ - ajyf + ia{^{a{^ + 2n){ai^ - 2N) = 0. (3.9) 

The boundary condition is aiv(t) ~ {uN /t) {1 + 0{t~^)) as t — > +oo. Note that 
Eq. (3.9), and therefore Eq. (3.8), contain the replica index n as a parameter. 

By derivation, Eq. (3.8) holds for n positive integers only and, strictly speaking, 
there is no a priori reason to expect it to stay valid away from n G Z"*". It can be shown, 
however, that it is legitimate to extend Eq. (3.8), as it stands, beyond n G Z+ and 
consider this extension as a proper analytic continuation to n G we are looking for 
(the reader is referred to [Kan02] for a detailed discussion) . 

As the result, the fermionic replica limit Eq. (1.13) can now safely be implemented 
to bring, via Eq. (2.6), the average density of eigenlevels [Meh60] 

^^'^ = 2^r(7V)0F [H'NmN-i{t) - Hr,ie)K_,{e)] (3.10) 

expressed in terms of Hermite polynomials. This is the famous GUE result firmly 
established by other methods [Meh04, Guh91]. Technically, the derivation of Eq. (3.10) 
is based on the small-n expansion of the Hamiltonian representation [Nou04] of the 
fourth Painleve transcendent. The details of this somewhat cumbersome calculation 

ft See also Chapter 9 of this Handbook. 

§ In the original paper [Kan02], Eq. (3.9) appears to have incorrect signs in front of n and N. I thank 
Nicholas Witte for bringing this fact to my attention. 
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can be found in [Osi09] where the nonperturbative result Eq. (3.10) is also re-derived 
within the bosonic variation of the replica trick. 

3.1.3. Brief summary 

The above treatment was largely based on a wealth of 'ready-for-use' results (Andreief- 
de Bruijn formula, Darboux theorem, and a connection between the Toda Lattice and 
Painleve transcendents) which surprisingly well fitted our goal of a nonperturbative 
evaluation of the particular replica partition function Eq. (3.1). Since existence of such 
an 'easy way' is clearly the exception rather than the rule, a regular yet flexible formal- 
ism is needed for a nonperturbative description of a general class of replica partition 
functions. 

3.2. The T function theory of replicas ((5 — 2) 

In this section, we outline such a regular formalism [Osi07, Osi09] tailor-made for an 
exact analysis of both fermionic and bosonic zero-dimensional replica field theories 
belonging to the broadly interpreted (5 — 2 Dyson symmetry class. 

3.2.1. From replica partition function to r function 

Let us concentrate on the fermionic and/or bosonic zero- dimensional replica field theo- 
ries whose partition functions admit the eigenvalue representation 

{^^n\<i))= [ nc^A,r(^;A,)e-^"(^^)A^(A), n e Z+. (3.11) 

Here Ki(A) is a 'confinement potential' which may depend on the replica index ±n; 
r(<;; A) is a function accommodating relevant physical parameters <; = (^i, • • •) of the 
theory (e.g., energies in the multi-point spectral correlation functions). In order to treat 
the fermionic and bosonic replicas on the same footing, the integration domain D was 
chosen to be f 

r 

'D^[j[c2j-uC2j]. (3.12) 
J=l 

To determine the replica partition function {Z^\q)) nonperturbatively, we adopt 
the 'deform-and-study' approach, a standard string theory method of revealing hidden 
structures. Its main idea consists of 'embedding' {Z^\q)) into a more general theory 
of T functions 

ri'H<;; t) = -,l i\d\e r(<j; A,) e-^"-^(^^) e"^*^^'^) A^(A) (3.13) 

t Notice that D = [—1,-1-1] for (compact) fermionic replicas, and D = [0,-|-oo) for (noncompact) 
bosonic replicas. A more general setting Eq. (3.12) does not complicate the theory. 
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which posses the infinite-dimensional parameter space t = (ti,t2,-'") arising as the 
result of the t-deformation 

oo 

v{t;X)^J2^,X^ (3.14) 

of the confinement potential. The auxiliary parameter s is assumed to be an integer, 
s e Z. Studying the evolution of r functions in the extended (n, s, t, <j) space allows 
us to identify the highly nontrivial, nonlinear differential hierarchical relations between 
them. Projection of these relations, taken at s = 0, onto the hyperplane t — 0, 

(2;W(.^))=n!r(^)(^;t) , (3.15) 

s=0,t=O 

will generate, among others, a closed nonlinear differential equation for the replica 
partition function {Z^n^\q)). Since this nonperturbative equation appears to contain 
the replica (or hierarchy) index n as a parameter, it is expected [Kan02] to serve as a 
proper starting point for building a consistent analytic continuation of {Z^\q)) away 
from n integers. 

Having formulated the crux of the method, let us turn to its exposition. The two key 
ingredients of the exact theory of r functions are (i) the bilinear identity [Dat83, Adl95] 
and (ii) the (linear) Virasoro constraints [Mir90]. 

3.2.2. Bilinear identity and integrable hierarchies 

The bilinear identity encodes an infinite set of hierarchically structured nonlinear dif- 
ferential equations in the variables {tj}. For the model introduced in Eq. (3.13), the 
bilinear identity reads [Adl95, Tu96, Osi07, Osi09]: 

-rt^-\t' - [^-^])^S^^^Si^ ) = 0. (3.16) 

Here, a e IR is a free parameter; the integration contour Goo encompasses the point 
2; = oo; the notation t ± [z~^] stands for the infinite set of parameters {tj ± /j}; for 
brevity, the physical parameters <j were dropped from the arguments of r functions. 

Being expanded in terms of t' — t and a, Eq. (3.16) generates various integrable 
hierarchies. One of them, the Kadomtsev-Petviashvili (KP) hierarchy in the Hirota 
form f 

1 AD,rW(t) orW(t) = s,^,{[D])Tj:\t)orj:\t) (3.17) 

I In Eq. (3.17), the j-th component of the infinite-dimensional vector [D] equals j~^Dj; the functions 
Sk {t) are the Schur polynomials [Mac98] defined by the expansion 



00 



exp ^tjX^ = ^x''sk{t). 
The operator symbol Dj f{t) o g(t) stands for the Hirota derivative f{t + x) g{t — x) \ 
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{k > 3) is of primary importance for the exact theory of rephcas. The first nontrivial 
member of the KP hierarchy reads 

S+^-^sSl) (|lo«-i"fc*))^-0. (3.18) 

In what foUows, it wiU be shown that its projection onto s = and t = [Eq. (3.15)] 
gives rise to a nonhnear differential equation for the replica partition function {Z^\<;)). 

3.2.3. Virasoro constraints 

Since we are interested in deriving a differential equation for {7}it\(;)) in terms of the 
derivatives over physical parameters {<;j}, we have to seek an additional block of the 
theory that would make a link between the {tj} derivatives in Eq. (3.18) taken at t = 
and the derivatives over physical parameters {<;j}. The study [Adl95] by Adler, Shiota, 
and van Moerbeke suggests that the missing block is the Virasoro constraints * which 
reflect the invariance of the r function [Eq. (3.13)] under a change of the integration 
variables. 

In the present context, it is useful to demand the invariance under an infinite set 
of transformations 

m 

A,^/., + e/.fV(/^,)n(/^.-4), ?>-l, (3-19) 

k=l 

labeled by integers q. Here, e > 0, the vector c' is d — {ci, • • • , C2r} \ {±00, i^} with i< 
denoting a set of zeros of /(A), and m = dim (c'). The function /(A) is, in turn, related 
to the confinement potential Vn-s{\) through the parameterisation 

^^fm' 9{>^) = Y.hX\ f{X) = Y.a,X' (3.20) 

■'^ ' fc=0 fc=0 

in which both g{X) and / (A) depend on n — s as do the coefficients and in the above 
expansions. The transformation Eq. (3.19) induces the Virasoro constraints [Osi07] 

£j'(t) + £j(<j;t)]r(^)(^;t) = 0, (3.21) 

where the differential operator 

= E E sra-ki-pjc')) aet,+,+eit) - be- ) , (3.22) 

e=o k=o ^ q+k+e+i/ 

acting in the t-space, is expressed in terms of the Virasoro operators § 

^.w = E.*.^+gs;s-. (3-) 

* Sec also Chapter 10 of this Handbook. 

§ Equation (3.23) assumes that d/dto is identified with the multiplicity of the matrix integral in 
Eq. (3.13), d/dto = n. 
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obeying the Virasoro algebra 

i^p, ^q] = {P- Q)^p+q, P,Q> -1- (3.24) 

The notation Sfc(— p^(c')) stands for the Schur polynomial and p„^(c') is an infinite 
dimensional vector 

pjc') = (^tr^(c'), ^tr^(c')', • • • , ^tr„(c')^ (3.25) 

where trm(c')'' = 

While very similar in spirit, the calculation of ^^{t), the second ingredient in 
Eq. (3.21), is more of an art since the function r(<j;A) in Eq. (3.13) may significantly 
vary from one replica model to the other. 

3.2.4- From r function to replica partition function 

Remarkably, for t = 0, the two equations [Eqs. (3.18) and (3.21)] can be solved jointly 
to bring a closed nonlinear differential equation for {Z^\q)). It is this equation which, 
being supplemented by appropriate boundary conditions, provides a truly nonperturba- 
tive description of the replica partition functions and facilitates performing the replica 
limit. 

3.3. Exact (bosonic) replicas at work: Density of eigenlevels in the chiral GUE 
3.3.1. Definitions 

To see the above formalism at work, let us consider the chiral Gaussian Unitary Ensem- 
ble (chGUE) of A/" X TV random matrices 



= ,,,t n (3-26) 




known to describe the low-energy sector of SU(Ai'c > 3) QCD in the fundamental 
representation ^ [Ver93]. Composed of rectangular x tir random matrices W with 
the Gaussian distributed complex valued entries 

iVE2 



/ 271 \"^"« 



-trW^W 



(3.27) 



where = ni + nn, the matrix !Kd has exactly u = Inj? — ulI zero eigenvalues identified 
with the topological charge u; the remaining eigenvalues occur in pairs {±Aj}; the 
parameter E denotes the chiral condensate. 

^ See also Chapter 32 of this Handbook. 
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3.3.2. Nonperturbative calculation of bosonic replica partition function 

To determine the (microscopic) spectral density from the bosonic rephcas, we define 
the rephca partition function 

{Z^-\,)) = {det-'^{<; + tn^))^ (3.28) 

(the angular brackets denote averaging with respect to the probability density Eq. (3.27)) 
and map it onto a bosonic field theory. In the half-plane D\c q > 0, the partition function 
{Z'h\q)) reduces to [DalOl, Fyo02] 

Zi-\u) = I dQr, det "-^Qn exp [-|Tr(0„ + Q-')] , (3.29) 

where the integration domain §„ spans all n x n positive definite Hermitean matrices 
Qn- Equation (3.29) was derived in the thermodynamic limit N ^ oo with the spectral 
parameter u; — qNT, being kept fixed (9lea; > 0). 

Spotting the invariance of the integrand in Eq. (3.29) under the unitary rotation 
of the matrix Q„, one readily realises that Zn~''(c<j) belongs to the class of r functions 
specified by Eq. (3.13) where D is set to 1R+, the potential Vns is 

Vn-s{X) ^{n-s-u) log X, (3.30) 

and r(^; A) is replaced with 



r(a;; A) = exp 



2 V A 



(3.31) 



This observation imphes that the associated r function Tn\uj; t) satisfies both the first 
KP equation (3.18) and the Virasoro constraints Eq. (3.21) with [Osi07] 

t^{t) ^t^^,{t) + {u-n + s)-^, (3.32) 

Projecting Eq. (3.18) taken at s = onto t = 0, and expressing the partial derivatives 
therein via the derivatives over uj with the help of Eqs. (3.21), (3.32) and (3.33), we 
conclude that 

hnico) = -^\ogZl-\u;) (3.34) 
obeys the differential equation [Osi07] 

.» + 5,._(4^i±i(^),U6(Ay^ 

CU \ UJ J 

+^^^^^^^^^/in - -^{Kf + -KK + ^ = (3.35) 

UJ^ LO^ UJ UJ^ 

that can be reduced to the Painleve III. 

Considered together with the boundary conditions hn{uj — > 0) ~ —nv/uj and 
hn{uj — > oo) ~ —n — n^/{2uj), following from Eq. (3.29), the nonhnear differential 
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equation Eq. (3.35) provides a nonperturbative characterisation of the bosonic rephca 
partition function * Zi, \oo) for all n E Z"*". 

3.3.3. Implementing the replica limit 

To pave the way for the replica calculation of the resolvent g{u!) determined by the 
replica limit 

Id 1 
g{uj) = -\im-—\ogZl^\uj) = -\im-hn{oj), (3.36) 

n^O n OUJ n^O n 

one has to analytically continue hn{uj) away from n integers. Previous studies [Kan02, 
Kan05] suggest that the sought analytic continuation is given by the very same Eq. (3.35) 
where the replica parameter n is let to explore the entire real axis. This leap makes the 
rest of the calculation straightforward. Representing hn{uj) in the vicinity of n = as 
hn{uj) — ^'^in^dpi^^), we conclude that g{u!) — — ai(a;) satisfies the equation 



g'" + 2a;V - (l + 4i^^ + 4a;^) uj^ + (1 - 41/^) ^ = 0. (3.37) 

Its solution, subject to the boundary conditions consistent with those specified below 
Eq. (3.35), 

g{u) ^- + +00 L(a;) K{u;) + h+i{uj) K_i{u;)] , (3.38) 
brings the microscopic spectral density g{u!) — 7r~^£He g{iu! + 0) in the form 

g{u) = u6{u) + I [j^iu) - J,.i{oo)J,+i{uj)j . (3.39) 

In the above formula (which is one of the celebrated RMT results originally obtained 
in [Ver93] within the orthogonal polynomial technique), the function J,y denotes the 
Bessel function of the first kind, whilst I,y and Ki^ are the modified Besscl function of 
the first and second kind, respectively. Let us stress that the approximate, saddle point 
approach to bosonic rephcas [DalOl] fails to produce Eq. (3.39). 

4. Concluding remarks 

Concluding this brief excursion into intcgrable theory of replica field theories, we 
wish to mention another important development due to Splittorff and Vcrbaarschot 
[Spl03, Spl04], not reviewed here for a lack of space. These authors showed that 
nonperturbative results for various RMT correlation functions at /3 — 2 can be derived 
by taking the replica limit of the graded Toda Lattice equation whose positive {n e Z+) 

* Interestingly, the fermionic replica partition function [Leu92, Smi95, DalOl] 

{w) = / DU„ det "Un exp [^Tr(U„ + Uj, 

admits a similar representation in terms of the Painleve III. Specifically, the function 
= {d/duj) \ogZ\a\'jj) can be shown [Kan02] to satisfy the very same Eq. (3.35) supplemented 
by the boundary conditions /„(a; u)q) = h-n{oJ coq), where cuq = {0, oo}. 
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and negative (n G Z~) branches describe fermionic and bosonic replica partition 
functions, respectively. Being a supersymmetric in nature, this approach greatly 
simplifies calculations of spectral correlation functions through a remarkable fermionic- 
bosonic factorisation. For further details, the reader is referred to the original papers 
[Spl03, Spl04], the lecture notes [Ver05], Chapter 32 of this Handbook, and the tutorial 
paper [Osi09]. 

Certainly, more effort is needed to accomplish integrable theory of zero-dimensional 
replica field theories. In particular, its extension to the /3 = 1 and (3 = 4 Dyson sym- 
metry classes is very much called for. 
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